WL-TP-92-015 


Progress  Towards  the  Development  of  Transient 
Ram  Accelerator  Simulation  as  Part  of  the  USAF 
Armament  Directorate  Research  Program 

Randall  P.  Drabczuk,  Captain,  USAF 

Wright  Laboratory,  Armament  Directorate 
Analysis  and  Strategic  Defense  Division 
Electromagnetic  Launcher  Technology  Branch 
101  W.  Eglin  Blvd.,Ste219 
Eglin  Air  Force  Base  FL  32542-6810 

G.  Rolader 

Science  Applications  inti.  Corp. 

Shalimar  FL  32547 

S.  Dash,  N.  Sinha,  B.  York 

Science  Applications  Inti.  Corp. 

Ft.  Washington  PA  19034-3211 


JANUARY  1993 


FINAL  PAPER  FOR  PERIOD  JANUARY  -  JUNE  1992 


or/. 

®*-®CTK 

FEB02I993 

6 


93-01854 


WHinniiic, 


Of>% 


Approved  for  public  release;  distribution  is  unlimited. 


WRIGHT  LABORATORY,  ARMAMENT  DIRECTORATE 

Air  Force  Materiel  Command  I  United  States  Air  Force  I  Eglin  Air  Force  Base 


98  2  1  070 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  in  connection  with  a  definitely  Government-related 
procurement,  the  United  States  Government  incurs  no  responsibility  or  any 
obligation  whatsoever.  The  fact  that  the  Government  may  have  formtdated  or 
in  any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be 
regarded  by  implication,  or  otherwise  as  in  any  manner  construed,  as  licensing 
the  holder,  or  any  other  person  or  corporation;  or  as  conveying  any  rights  or 
permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  in  any 
way  be  related  thereto. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 

The  Public  Affairs  Office  has  reviewed  this  report,  and  it  is  releasable  to  the 
National  Technical  Information  Service  (NTIS),  where  it  will  be  available  to 
the  general  public,  including  foreign  nationals. 

FOR  THE  COMMANDER 


Even  though  this  report  may  contain  special  release  rights  held  by  the 
controlling  office,  please  do  not  request  copies  from  the  Wright  Laboratory, 
Armament  Directorate.  If  you  qusJify  as  a  recipient,  release  approval  will  be 
obtained  from  the  originating  activity  by  DTIC.  Address  your  request  for 
additional  copies  to: 


Defense  Technical  Information  Center 
Cameron  Station 
Alexandria  VA  22304-6145 


If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing 
list,  or  if  your  organization  no  longer  employs  the  addressee,  please  notify 
WL/MNSH,  Eglin  AFB  FL  32542-6810,  to  help  us  maintain  a  current  mailing 
list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notice 
on  a  specific  document  requires  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  ii  estimated  to  average  i  hour  per  response,  including  the  time  for  reviewing  instruaions.  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  arm  reviewing  the  coileaion  of  information.  Ser>d  comments  reaarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  I21S  Jefferson 
Oavts  Highway,  Suite  1204.  Arlington,  VA  22202'4302.  and  to  the  O^iceof  Management  arm  Budget.  Paperwork  Reduction  Project  (0704^188),  Washington,  DC  20S0). 


1.  A6€NCY  USE  ONLY  (leave  biank) 


2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

January  1993  Final  January  to  June  1992 


4.  TITLE  AND  SUBTITLE  5.  FUNDING  NUMBERS 

Progress  Towards  the  Development  of  Transient  Ram 

Accelerator  Simulation  as  Part  of  the  USAF  Armament  PEi  62602F 

Directorate  Research  Program  PRt  2502 

6.  AUTHOR(S) 

Randall  P.  Drabczuk,  Capt.  USAF  VUt  04 

G.  Rolader.  SAIC 

S.  Dash,  N.  Sinha,  B.  York,  SAIC 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADORESS(ES) 


B.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


Same  as  Block  9  together  with 

Science  Applications  Inti.  Corp.,  Shalimar  FL  32579 
Science  Applications  Inti.  Corp.,  Fort  Washington  PA 

19034-3211 


WL-TP-92-015 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 

Wright  Laboratory,  Armament  Directorate 
Analysis  and  Strategic  Defense  Division 
Electromagnetic  Launcher  Technology  Branch  (WL/MNSH) 
101  W.  Eglin  Blvd.,  Ste  219 
EgUn  AFB  FL  32542-6810 


11.  SUPPLEMENTARY  NOTES 

Approved  by  PA  for  unlimited  release  in  June  1992.  Published  in  the 
AIAA/SAE/ASME/ASEE  28th  Joint  Propulsion  Conference  proceedings. 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

This  paper  describes  the  development  of  an  advanced  computational  fluid  dynamic 
(CFD)  simulation  capability  in  support  of  the  USAF  Armament  Directorate  ram 
accelerator  research  initiative.  The  state-of-the-art  CRAFT  computer  code  has  been 
specialized  for  high  fidelity,  transient  ram  accelerator  simulations  via  inclusion  of 
generalized  dynamic  gridding,  solution  adaptive  grid  clustering,  and  high  pressure 
thermochemistry.  Selected  ram  accelerator  simulations  are  presented  that  serve  to 
exhibit  the  CRAFT  code  capabilities  and  identify  some  of  the  principle 
research/design  issues. 


14.  SUBJEa  TERMS 

Ram  Accelerator,  Hypervelocity  Launcher,  Oblique  Detonation 
Wave,  Chemical  Propulsion 


IS.  NUMBER  OF  PAGES 
20 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  118.  SECURITY  CLASSIFICATION  119.  SECURITY  CLASSIFICATION  I  20.  LIMITATION  OF  ABSTRACT 

OF  REPORT  I  OF  THIS  PAGE  I  OF  ABSTRAO  I 


UNCLASSIFIED 


NSN  7540-01-280-5500 


UNCLASSIFIED 


UNCLASSIFIED 


Standard  Form  298  (Rev  2*89) 
Pr^ribed  by  ansi  Std  2)9-18 
298-102 


PREFACE 


This  work  was  funded  by  WL/MNSH  of  the  Armawent  Directorate  at  Eglin  AFB 
FL  under  the  Ran  Accelerator  Research  progran.  Capt  Randall  P.  Drabczuk. 
WL/MNSH,  and  personnel  from  Science  Applications  International  Corporation 
(SAIC)  in  Shaliaar  FL  and  Ft.  Washington  PA  performed  the  work  during  the 
period  of  January  1992  to  June  1992. 

This  technical  paper  was  presented  at  the  AIAA/SAE/ASME/ASEE  28th  Joint 
Propulsion  Conference  in  Nashville  TN  on  6-8  July  1992. 


iii/iv  (Blank) 


Wowasslon  For 


WIS  GRAftI 
DTIC  tab 
Unanrouiiped 
Just  LT  lest  Ion _ _ 


□ 

□ 


By— _ _ 

Dlgtrlbutlon/ 


fDlbt 


Availability  Codes 
j2Vall  and/or 


Spaoial 


DTIC  QUALITY  INSPECTED  3 


PROGRESS  TOWARDS  THE  DEVELOPMENT  OF  TRANSIENT  RAM  ACCELERATOR  SIMULATION 
AS  PART  OF  THE  U.S.  AIR  FORCE  ARMAMENT  DIRECTORATE  RESEARCH  PROGRAMf 

N.  Sinha  *  B.  J.  York.  **  and  S.  M.  Dash*** 

Science  Applications  International  Corporation 
Propulsion  Fluid  Dynanucs  Division 
Fort  Washington,  Pennsylvania 

R.  Drabczuk**** 

Wright  Laboratory 
Araamenl  Directorate 
Eglin  AFB,  Florida 

G.  E.  Rolader***** 

Science  Af^lications  International  Corporation 
Shalimar.  Florida 


APsma: 

This  paper  describes  the  development  of  an  advanced 
conoputational  fluid  dynamic  (CFD)  simulation  capability  in 
suiqport  of  the  U.S.  Air  Force  Armament  Directorate’s  ram 
accelerator  research  initiative.  The  state-of-the-art  CRAFT 
oon^niter  code  has  been  specialized  for  high  fidelity,  transient 
ram  accelerator  stmulstions  via  inclusion  of  generalize  dyiuunic 
gridding,  solution  adaptive  grid  clustering,  high  pressure  drermo- 
chemistry,  etc.  Selected  ram  accelerator  simulations  are  pre¬ 
sented  whidi  serve  to  exhibit  the  CRAFT  code's  capabilities  and 
identify  some  of  die  principal  research/design  issues. 

INTRODUCTION 

At  an  earlier  AIAA  meeting,*  a  new  propulsive-ouented 
Navier-Stokes  ^S)  code,  CRAFT,  was  introduc^  aduch  was  an 
extendon  of  ^  TUFF  aerodynamic  NS  code  of  Molvik  and 
Meifcle.^  Applicatirms  of  CRAFT  to  various  combusting, 
multiphase  flw  problems  were  exhibited.  In  this  pqier,  the 
siaiulation  oi  lam  accdeiator  flowfields  widi  CRAFT  wfll  be 
discussed  -  a  role  undertaken  in  support  of  new  research  initia¬ 
tives  at  the  Armament  Directorate  of  Wright  Labs,  U.S.  Air 
Force. 

The  Aimament  Directorate  Wright  Labcuatoiies 
(WL/MN),  located  at  Eglin  Air  Force  Base,  in  conjunction  with 
dm  Air  Force  Office  of  Scientific  Reseaidi  (AFOSR),  has 
initiated  a  multi-year  research  program  to  advance  the  under¬ 
standing  of  the  basic  physics  of  the  ram  accdeiator.  The 
researA  effort  has  three  principal  areas  of  focus: 
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(1)  dengn  and  construction  of  a  large  bore  ram  accelerator 
test  ftcility  at  Eglin  AFB; 

(2)  investigation  of  diagixistic  tedmiqnes;  and, 

(3)  devdopment  and  execution  of  computational  fluid 
dynamic  (CFD)  models. 

The  design  of  the  ram  accdeiator  researdi  fiidlily  to  be 
constructed  at  Eglin  AFB  has  been  initiated.  Present  plans  call 
for  the  qrstem  to  include  a  conventional  powder  gun  injector,  an 
evacuated  vent  charnber,  an  dghteen  meter  long  ram  accdeiator 
section  with  a  90  nun  bore  diameter,  and  a  catdi-tank. 

The  diagnostic  research  will  focan  on  schlieren,  x-rqr 
diadowgiapb,  emission  ^ectroicopy,  and  a  single-poise  imagms 
(Rayleigb  scattering).  The  seUieren  todoiiqoe  depends  on  fte 
diange  in  Hhumnation  due  to  dm  dtaagea  in  lefnctive  index 
(densiQr  grarSents),  wdiicii  are  made  visible  by  mtroduemg 
subadiary  dements  into  dm  optical  qrstem,  suA  as  a  knife  e^e, 
resulting  in  the  ippearance  of  dm  intenafy  variations  ia  dm 
location  of  dm  focussed  image.  T1mx-ny  tednnquadqieodtoo 
dm  diange  in  dm  iHununatioo  due  to  absiMptifla  akatg  dte  line  of 
dght  by  dm  gas  constituents,  in  a  — w"**  amQar  to  optied 
diadowgnphy,^^  wfaidi  examines  fee  intensify  variations  in  a 
plaimwhi^  is  not  the  plaim  of  the  focussed  image.  Tbeinlerest 
in  nang  x-rays  deiivcB  from  a  need  to  have  a  nmdiod  feat 
operates  successfully  at  hi^  pressuiea.  X-ray  ahadowgrfybiy 
htt  been  utilized  to  observe  denafy  gradfents  in  dm  flow  as  dm 
pressure  increases  up  to  values  ia  dm  kilobar  lagime.  Eaasaon 
spectroscopy  will  be  enployed  to  infer  flowfidd  temperatare  ad 
pressure,  as  weD  as  conapositian.  An  antidpated  difficulfywife 
emissioo  fyectra  is  that  the  opadfy  of  dm  flow  fystem  mqr 
increase  rapidly  with  iacreasiag  pressure.  The  lesearfe  wiD 
idatify  dm  pressure  rfyime  over  whidi  diis  lediaique  will  be 
useful.  The  final  diagnostic  techaique  to  be  invmtigated  is 
dagle-pulse  planar  imaging.  The  primary  focus  wiO  be  ca 
Rayld^  scatteriag,  but  potatially  both  Ranam  and  Tbompsa 
scattering  will  be  utilized. 

High  fiddify  CFD  sianilations  win  be  a  initial  part  of  dm 
overdl  ram  accdeiator  researdi  program.  CFD  rimulatkmswiB 
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be  used  for  conducting  tradeoff  studies  to  support  the  in-house 
lest  facility  design  end  develt^ment.  A  second  objective  of  the 
cosDputationsl  fluid  dynsnucs  task  is  to  stfport  the  development 
of  the  diagnostic  tec^iques  and  data  inteq)reUtion.  The  flnal 
and  crucial  longer  term  objective  is  to  develop  a  detailed  under¬ 
standing  of  the  complex  fdiysical  and  theimocbemical  processes 
governing  the  ram  accelerator  concept.  Ibis  paper  will  provide 
a  background  of  the  CFD  technique  ^ing  utilized  to  support  this 
research  initiative  and  present  preliminary  flndings. 

Earlier  simulations  of  ram  accelerator  flowflelds  have 
primarily  emphasized  a  ‘snapshot*  prediction  of  the  ram  flow- 
field,  e.g.,  a  steady-state  simulation  at  a  fixed  trajectory  point. 
Such  studies  include  the  work  of  Nusca  at  BRL,*  Yungster  and 
Eberhardt  at  the  University  of  Washiqgton,***^  Yungster  at 
NASA  Lewis,*  Imlay  ct  al.  at  AMTEC,*^  Humphry  and 
Sobota  at  APRI,*^  and  Li,  Kailasanath  and  Oran  at  NRL.*^  The 
work  of  Nusca,  Yungster  et  al.,  and  Imlay  et  al.  has  involved 
the  use  of  modem  fim'te-vohune  upwindfimplicit  numerics  (e.g., 
higher-order  Roe/TVD).  Nusca  has  used  uncoupled,  global 
chemistry  in  an  overlaid  numner.  Yungster  and  Imlay  et  al. 
have  used  loosely-coupled  (diagonalized)  finite-rate  chemistry. 
Li  et  al.  have  utilized  explicit  FCT  numerics  with  time-qilit 
global  chemistry.  With  regard  to  turbulence  modeling,  Nusca, 
Yungster  and  Iiday  et  al.  have  employed  algebraic  models  (e.g., 
Baldwin-Lomax),  vdiile  Li  et  al.  have  not  utilized  a  turbulence 
model  (their  work  is  Euler  or  viscous  with  ‘boosted*  laminar 
viscosity). 

A  workshop  was  recently  held  at  AFOSR  wifli  a  primary 
purpose  of  identifying  principal  research  issues  and  requirements 
for  high  fidelity  CFD  simul^on  of  ram  accelerator  flowfields. 
The  workshop  group  felt  fast  a  full  transient  simulatioo  was 
requisite  for  studying  the  basic  physics  and  thermochemistry  of 
lam  accderators.  Very  little  remardi  has  been  done  in  tftis 
area.  The  only  work  to  date  using  staleH)f-lhe-art  upwind/im¬ 
plicit  numerics  to  analyze  transient  hypervelocity  interior  ballistic 
flows  is  that  of  Sinha,  Dadi  el  al.  as  applied  to  Electro-Tbermo- 
Cbemical  (ETC)  guns.**  Alternate  transient  interior  ballistic 
flow  simulalion  work  for  ram  acceleiatais,  for  ETC  guns,  and 
for  liquid  propellaat  guns,  has  primarily  invdved  die  use  of 
structured  and  unstructured  explicit  methods  (vriiich  caimot 
efficieady  treat  thin  boundary  layers).  The  principal  researdi 
areas  identified  at  the  workshop  included  numerics  (e.g., 
dynamic  adaptive  gridding  and  advanced  impIicilAqiwind 
algorithms),  turbulence  modeling  in  a  transieat  environment 
(e.g.,  LES  simulation),  3D  effects  (e.g.,  vortex  diedding  from 
fins)  and  combustion  phenomenology  in  the  high  Re,  high 
pressure  ram  acceleralor  environment 

The  methodology  in  CRAFT  is  quite  new  and  reflects  ‘state- 
of-the-art*  capabilities  in  this  country.  Principal  milestones  in 
the  evolution  of  CRAFT  are  summarized  m  Table  I.  Our  work 
with  earlier  NS  codes  such  as  PARCH*^’**  has  indicated  the 
requirement  to  develop  qiedalized  versians  for  specific  applica¬ 
tions.**'^  This  philosr^y  has  been  adhered  to  in  our  work 
with  CRAFT.  Specialize  research  versions  of  CRAFT  devel¬ 
oped  to  date  are  described  in  Table  II.  This  paper  will  review 
Ite  basic  methodology  available  in  CRAFT  (equations  and 


numerics)  supplementing  descriptions  provided  in  earlier 
publications.*'^  The  problem  qiecific  features  of  the  RAM 
version  will  be  emphasized. 

GAS-PHASE  EQUATIONS 

The  Reynolds- Averaged  Navier-Stokes  equations  integrated 
by  the  RAM  version  of  CRAFT  are  listed  in  Table  m.  They 
include  chemical  species  equations  for  the  qiecies  mass  fnctioos 
p;,  and  turbulence  model  equations  for  the  turbulent  kinetic 
energy,  k,  and  turbulence  dissipatiao  rate  parameter,  c.  The 
source  terms  in  D  provide  for  chemical  non-equilibrium  (6^, 
and  turbulence  non-equilibrium  (£-pc)  where  2  i*  the  production 
term  for  turbulent  kinetic  energy.  The  viscous  fluxes  indude 
full  cross-flow  stress  terms  wfaidi  is  requisite  fiv  solving  3D 
problems,  e.g.,  projectile  with  fins.  A  two-equation  kc  turbu¬ 
lence  modeH*  has  been  incorporated  in  a  strongly  coupM, 
implicit  manner.  Near-wall  effects  have  been  dealt  with  via  the 
low  Reynolds  number  formulation  of  Chiei^  and  the  recent 
modd  of  Rodi.^  High-speed  compressibility  issues  have  been 
accounted  for  via  the  upgrades  of  Dadi  et  al.^*  based  on 
‘compressible  dissipation*  concqits  introduced  by  Sarfcar^  and 
Zeman.^  The  CRAFT  code  «»"**«"«  generalized,  multi- 
component/multi-st^,  firute-rate  diemical  kinetics.  The  ^ledes 
production/destruction  rates  ate  determined  via  the  law  of  mass 
action  and  chemical  reaction  rate  constants  are  expressed  in 
Arrhenius  form.  The  code  accesses  a  JANNAF-baaed 
thermocbemical/tranqxMt  databank  whose  features  mimic  tiiose 
in  standardized  codes  such  as  SPF.^  The  Aitab«nV  contains 
thermodynamic  and  tranqiort  data  curve  fits  for  varied  qpedcs 
and  rate  constants  for  pertinent  diemical  reactions. 

NUMERICAL  PROCEDURES 

The  CRAFT  code  is  an  extended  propulsive-orienled 
derivative  of  the  TUFF  aerodynannc  code  vdiidi  was  formulated 
as  a  thin-layer  NS  code  within  a  finite-vcdume  fiamewoik,  with 
higher-order  upwind  invisdd  fluxes  based  on  the  approxinnte 
Riemann  solver  of  Roe.^  The  numerics  mcoipmsle  TVD 
procedures  outlined  by  Chdoavstdiy^  to  dlmlnate  spurious 
oscfllations  generally  assodsted  widi  higher  order  upwind 
schemes.  The  fluid  ^  qiedes  equations  ate  strongly  coupled, 
and  the  solution  is  made  fully  implidt  by  linearizing  aD  convec- 
tive/difluave  fluxes  as  wdl  as  the  chemical  source  term  (wife 
reflect  to  all  dqiendent  variaUes,  not  just  chemical  qpedes  as 
in  other  codes  of  fliis  type).  An  ADI  qdittiiig  is  cm|doyed.  A 
modified  Newton  iteration  is  employed  to  rfimiwte  linearization 
errors  in  the  iterative  process.  EUmiiudion  of  linearization/fac- 
torization  errors  was  recently  diowa  by  Ridder  and  Beddini***  to 
be  essential  for  accurately  rqiresenting  corriplex  wave  propaga¬ 
tion  pbenomeru  in  a  non-difiusive  marmer,  as  requisite  fat 
acoustic  applications  and  instaNliqr  phenomena. 

The  incorporation  of  nonequilibrinm  kinetics  has  followed 
the  finite-volume,  fully-inqilicit  philosophy  of  die  gas  dynamics. 
In  dqiarture  firom  earlier  generation  explicit  and  loosely  coupled 
schemes,  the  species  and  fluid  dynamic  equations  are  strongly 
coupled  and  simultaneously  solved,  readily  allowing  infbtmation 
transfer  between  the  two  sets  of  equations.  In  ccmtrast,  the 
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Table  I.  Evolution  of  Methodolosy  in  CRAFT  Code 


•  UPWIND/IMPLICIT  FINITE-VOLUME  METHODOLOGY  (ROE/TVD)  APPLIED  TO  EXTERNAL  AERODY¬ 
NAMIC  ^OWFIELDS  ^  1985-1988.  Reference:  Walters  and  Thomas,  ‘Advances  in  Upwmd  Relaxation  Methods,” 
State  of  the  An  Surveys  on  Computational  Fluid  Mechanics,  ASME  Pub.,  1988. 

•  EXTENSIONS  OF  ROE/TVD  METHODOLOGY  TO  MULTI-COMPONENT/REAL  GAS  FLOWS  WITH  CHEMI¬ 
CAL/THERMAL  NONEQUILIBRIUM  ~  1988-1989.  Reference:  Liu  and  Vinokur,  ‘Upwind  Algorithms  for  Genetal 
Dienno-Cheiiucal  Nonequilibrium Flows,”  AIAA  Paper  89-0201,  Jan.  1989. 

•  EXPLORATORY  STUDIES  OF  ROE/TVD  NUMERICS  FOR  UNSTEADY  FLOWS,  ACCURATE  W'AVE 
PROPAGATION  IN  DUCT  (RESONANCE  TUBE)  ~  1989.  Reference:  Ridder  and  Beddini,  ‘Time-Accunte  Finite- 
Volume  Method  for  Propulsion  Chamber  Flows,”  AIAA  Paper  89-25S4,  July  .1989. 

•  EXTENSIONS  OF  ROE/TVD  METHODOLOGY  FOR  PROPULSIVE  FLOWS  WITH  COMBUSTION  CHEMIS¬ 
TRY,  ADVANCED  TURBULENCE  MODELS  ~  1989-1991.  Reference:  Dash.  'Advanced  Computational  Models  for 
Analyzing  High-Speed  Propulsive  Flowfields,”  I990JANNAF  Propulsion  Meeting,  Oct.  1990. 

•  EXTENSIONS  OF  ROE/TVD  METHODOLOGY  FOR  UNSTEADY  ROCKET  PLUME  PROBLEMS  WITH 
GAS/PARTICLE  NONEQUILIBRIUM  --  1989-1991.  Reference:  Sinha,  Hosangadi  and  Dash.  ”T1ie  CRAFT  NS  Code 
and  Preliminary  Applications  to  Unsteady,  Reacting,  Multi-Phase  Plume  Flowfield  Problems,”  J992  JANNAF  Plume 
Meeting,  May  1991. 

•  EXPLORATORY  STUDIES  OF  ROEfTVD  METHODOLOGY  FOR  GUN  FLOWTIELDS  WITH  DYNAhOC  GRID 
EXTENSIONS  —  1990-1991.  Reference:  Hosangadi  and  Sinha,  ‘Development  of  a  Dynamic  Grid  Navier-Stokes  Solver 
for  Computing  Electro-Thermo-Chemical  Gun  Configurations,*  SAIC/TR-88,  Dec.  1990. 

•  EXTENSIONS  OF  ROE/TVD  METHODOLOGY  FOR  GAS/LIQUID  hflXTURES  FOR  ETC  GUN  FLOWFIELD 
SIMULATION  —  1992.  Reference:  Sinha,  Hosangadi,  York  and  Dash,  ‘First  Principles  Two-Dimensional  Modeling  of 
ETC  Interior  Ballistics,”  29th  JANNAF  Combustion  Meeting,  October  19^. 


Table  U.  Research  Versions  of  CRAFT  NS  Code 


CODE  NAME 

CRAFT/JR 

CRAFT/ETC 

CRAFT/LU 

CRAFT/RAM 

APPUCATION 

High.Spccd  Icl 

Rcfcatch 

ETC  Cun  Fkwficldi 

Rocket  Nozzlea/ 

Phimea 

Numerical  Reaearcb 

Ram  ArreWtaeor 
Flowfielda 

SPONSOR 

NASAAaRC,  HSRP 

BRL 

MiCOM 

iMcmal  Rcaeaicb 

AFOSR* 

WL/MNSH 

ESSIKSiH, 

ID/2DrAXV30 

1D/3D/AX1/30 

ltV3D/AXV3D 

IDGO/AXIOD 

IDOD/AXiaO 

CHEMISTRY 

Perfect  Get,  Two* 

Streeiii 

Imperfect  Oat,  Conibiia- 
tion  Chcmifiry,  Muki- 
Phaac  Chcfidairy 

SPF  Finilc-Ralc 
Cbcfliiaiiy.  Two- 
Stream,  P^eci  Gaa 

PertcclGaa 

Rjm  D^iled  Kmetica 

TURBUUENCE 

ke,  Comprenibility  Ex- 
tentkNit,  ARS  Exi^ 
tiont 

ke 

kc/Chien,  Comprata- 
RiSiy  ExtcMioae 

ke 

ke/CMaa,  ABvMcad 

racUons 

PARTICULATES 

None 

Equilibrated  Mixnife, 
Nonequilibiwm  Liquid 
PropellaiXa 

Noiwquilibriuro 

None 

None 

SOLUTION 

Implicil/Upwiad 
CRoe/TVD)  Strongly- 
Coupled  Fluid/Speciet/ 
Turbulence  (1  x  t 

Block  InveniorO, 
Time-Accurate  or 
Time-Atympioiie 

ImpCcittUpwind 
GtoerrVD)  on  Oynamie 
Grid,  Timc-Accuralc 
Solution,  Strongly-Cou¬ 
pled  Equationa,  Vari> 
able  Matrix  Size 

ImpBcilAJpa'iad 
(Roe/TVD)  on  Fixed 
Grid,  Time-Accurate 
or  Time-Aaympioiie, 
Strongly-Coupled 
Equationa,  Variable 
Matrix  Size 

LU  Upgrade  far 
RobuUMtt,  Fealer 
COoveigcnce 

CCFL  -  25-Slb 

Impfieil/Upwind 
<Roe/TVD}  on  Qy- 
namic  Grid,  Time-Ae- 
curate  Solutiaa, 
Stroogly-Coeplad 
EquaGona,  Variabk 
Matrix  Sixa 

NEW  WORK 

Exploriitg  New  Tuibii* 
lence  Modelt,  New 
Non-Reflcctive  BC 

Inclution  of  liquid 

Ftwae,  Gaa/Iiquid 
Inurfice,  Droplet  For- 
mttiM  Model 

Adaptive  Gridding  for 
Unaieady  Muhi-Tbaae 
Fkiwa 

Rewrite  of  Code 
Stiuemie  to  OptimiM 

LU  Storage  fcr  3D 

Adaptive.  Dynn^ 
Otid«i«  fa  OHUady 
FIowb,  LES 
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Table  m.  CRAFT  NS  Equations 


♦  /M  w  -  *  /M*«i  - 

-  ^  //(*..»  -  *  //Pi.ui  -  *i.uO<WC 


■  X(u,^v,»w^-»2|ni, 

P*  -  -KT,*ut.*VT^^wt^ 

P.  •  -CT,*ut^*vt^^wT. 
X  •  -%|» 


IN  THE  ABOVE  EQUATIONS.  U*.rA,  REPRESENT  THE  CURVCJNEAR  COORMNAIES  IN  THE  TRANSFORMED 

COMPUTATIONAL  SPACE.  THE  METRICS  f,  ARE  THE  COMPONENTS  OF  THE  CELL  FACE  NORMAL  7.  POtNITNO  IN  THE 

POSITIVE  tORECnONWnH  LENGTH  EQUAL  TO  CELL  FACE  AREA.  SMUARLY  THE  OIHBt  METRIC  QUANTITIES  ARE  THE 
COMPONEKTSOFTHEVECTORSmANDB.WHICHARETHENORMALVECTORSOFTHECORRESPONI>INOH,ANDa.  IHEFIMD 
PRWERTIES  ARE  GIVEN  BY  THE  DENSITY  0,  CARTESIAN  VELOCITY  COMPONENTS  (m  y.  w).  AND  THE  TOTAL  ENBtOY  PBt  UNIT 
VOLUME.*.  AND  SPECIES  CONCENTRAT1ONS,ai.THETRANSP0RT  PROPERTIES  AREOIVENBYiiBElNOTHELAMmARVlSOOSITY. 
V  BEINO  THE  RATIO  OF  SPECIFIC  HEATS.  AND  KBEINO  THE  THERMAL  CONDUCTIVITY.  THE  TURBULENCE  EQUATIONS  ARB 
DERIVEDD4TERMSOFTHETURBULENTKINEnCENEROY.k.ANDTHEDISSlPAT10NRATE,c.  THE OHCIRAVARIANT VELOCITY 
COMPONENTS  ARE  DEFINED  AS  FOLLOWS: 

0  -  ♦  ^V  ♦ 

V  ■  SB^V  ♦  BlyV  ♦  BIjW 

W  ■  n,®  ♦  lyr  ♦  i^w 


FOR  TURBULENT  FLOW  SIMULATION.  THE  LAMINAR  VISCOSITY  AND  THE  THERMAL  CONDUCTIVITY  ARE  REPLACED  BY  THE 
RESPECTIVE  SUMS  OF  LAMINAR  AND  TURBULENT  VALUES.  THE  k-*  TURBULENCE  MODEL  CONSTANTS  WHICH  RAVE  BEEN  U»D 
AREC  -  .09.C,  -  1.44 uNC,-  1.92. 


■  earlier  techniques  sepantely  solved  (he  tu’o  sets  of  equations 
which  reduces  the  computational  effort  but  severely  inhibits  the 
flow  of  information.  This  can  be  a  severe  handicap  for  strongly 
combusting  flows,  e.g.,  blast  waves,  gun  chamber  flowfields, 
etc.,  where  signiflcant  coupling  exists  between  the  chemical 
eelease  of  energy  and  flowfield  pressurization.  A  major 
breakthrough  has  been  the  pioneering  work  of  Vinokur  and  Liu^* 
in  reformulating  the  original  'perfect  gas*  Roe  scheme  and  the 
.associated  Reinuim  problem  to  account  for  multi-component, 
real  gas  flow  with  chemical  and  thermal  nonequilibrium. 

The  extension  of  CRAFT  to  include  a  two^uation  kc 
turbulence  model  is  a  straightforward  extension  of  the  basic 
finite-volume  based  upwtnd/TVD  methodology.  The  turbulence 
model  equations  are  cast  as  a  part  of  the  strongly-coupled 
equation  set  for  enhanced  accuracy  and  robustness.  The  implicit 
treatment  is  extended  to  the  turbulence  nonequilibrium  source 
terms  to  eliminate  'stifhiess.'  This  is  requisite  for  handling  the 
low  Reynolds  number  near  wall  terms.  The  linearization  of 
source  terms  is  a  somewhat  inexact  procedure  and  several 
strategies  have  been  developed.  Instead  of  a  brute  force 
mathematical  treatment,  the  linearization  is  being  guided  by 
turbulence  behavior  typical  of  low  Reynolds  number  near  wall 
regions  and  high  Reynolds  number  free  shear  zones.  Various 
methods  are  bemg  explored  to  identify  the  optimal  technique  for 
propulsive/plume  flowfields.  In  addition,  the  selection  of 
qiprt^riate  background  levels  for  freestream  fluctuation  intensity 
and  length  scale  have  been  found  to  be  sensitive  parameters  and 
must  be  judiciously  specified. 

DYNAMIC  GRIP  UPGRADES 

As  a  preliminary  step  towards  developing  a  gun  capability. 
CRAFT  was  upgraded  to  allow  for  unsteady  simulations  with  a 
moving  projectie  and  thus  an  ever  expanding  and/or  moving  gun 
chamber  computational  domain.  The  generalized  dynamic  grid 
formulation  incorporated  enhances  (he  capabilities  of  the  code  to 
«itniil«tA  a  wide  variety  of  other  unsteady  phenomena,  e.g.,  the 
unsteady  interactions  between  deformable  or  non-slationary 
boundaries  and  fluid  phenomena  which  result  in  time-varying 
flowfields.  Dynamic  grids  are  also  applicable  to  problems  which 
require  a  grid  that  continuously  adapts  to  moving  discontinuities. 
This  feature  is  requisite  in  transient  problems  where  strong 
shocks  can  oonvect  rapidly  through  the  flowfield. 

The  dynamic  grid  capability  has  been  formulated  in  a 
manner  consistent  with  the  finite-volume  philosophy  of  the 
CRAFT  code.  We  have  closely  followed  the  generalized 
procedure  outlined  by  Vinokur.^^  Table  IV  summarizes  this 
methodology,  while  additional  details  are  given  in  Ref.  33.  To 
briefly  describe  the  procedure,  the  unsteady  movement  of  the 
physical  grid  is  treated  purely  as  a  geometric  quantity.  There¬ 
fore,  unlike  a  finite-difference  formulation  where  (he  contravar- 
iant  velocity  is  altered  such  that  the  computational  grid  remains 
invariant  in  qtace  and  time,  the  computational  grid  in  a  fi¬ 
nite-volume  formulation  varies  with  time.  The  additional  flux 
due  to  the  motion  of  the  cell  face  is  included  in  the  flux  evaula- 
tioo  in  an  upwind  fashion.  The  conservation  equations  incorpo¬ 
rate  die  chuge  in  the  flux  resulting  from  the  time-varying  cell 


volume  as  a  source  term.  The  variation  of  the  grid  with  time 
has  to  be  specified  depending  on  the  application.  For  a  moving 
boimdary,  the  movement  of  the  boundary  has  to  be  related  to  the 
force  balance  on  it.  On  the  other  hand,  for  a  grid  adapting  to 
convecting  discontinuities,  the  criteria  for  the  adaption  (e.g. 
pressure  or  density  gradient)  must  be  specified. 

RAM  ACCELERATOR  STUDIES 
Steady-State  Simulations 

It  was  imperative  to  first  establish  a  baseline  capability  at 
the  steady-sute  axisymmetric  level  to  establish  that  CRAFT 
provided  jpredictions  in  accord  with  those  of  other  investiga¬ 
tors.**' With  this  objective  in  mind,  predictions  were 
generated  for  a  projectile  that  was  recently  tested  at  the  Univer¬ 
sity  of  Washington’s  38  nun  ram  accelerator  test  facility. 

The  projectile  dimensicms  are  shown  in  Figure  1.  The  tube 
was  pressurized  to  39  atm  with  a  premixed  combustible  gas 
mixture  (2.75  CH4  -f  5.8  Nj  *<'  2  O^.  CRAFT  predictions 
were  obtained,  at  ^e  Euler  non-reacting  limit,  for  a  trajectory 
point  corresponding  to  a  projectile  flight  Mach  number  of  3.6. 
A  301  X  71  curvilinear  body  fitted  grid  (shown  in  Figure  2)  was 
utilized.  Figures  3  and  4  show  ccnqruted  pressure  and  Madt 
number  contours,  respectively.  The  overall  flow  structure  is 
well  rqrresented  and  in  accord  with  predictions  obtained  by  the 
previous  investigators  utilizing  conqtarable  numerics.  The 
oblique  shock  structure  is  initiated  at  the  projectile  nose  cone  and 
causes  ram  compression  of  die  reactants  via  a  sequence  of 
oblique  shocks.  The  shock  system  weakens  as  it  interacts  with 
the  expansion  fan  emarutmg  from  the  projectile  afterbody 
shoulder.  Behind  the  projectile,  a  bigb-speed  wake  structure  is 
obtained,  characterized  by  a  wake  leattachment  diode. 

High-IVessure  Compressibnity  Effects 

Ram  accelerator  flowfields  ate  characterized  by  extremdy 
high  pressures,  caused  by  rapid  combustion  of  pie-pressurized 
semi-stoichiometric  fiiel-air  mixtures  under  near-detonative 
dtuations.  Pressures  in  the  nd^borbood  of  1000  atmoqiheres 
are  not  unusual.  Under  such  extreme  pressure  oonditioos,  high 
levels  of  compressibility  are  prevalent  and  must  be  accounted  for 
during  analysis. 

The  baseline  version  of  CRAFT  was  designed  for  a  multi¬ 
component,  thermally  perfect  (but  calorically  inqierfed)  gaseous 
mixture  and  the  perfect  gas  equation  of  state  was  assumed,  h 
the  RAM  version,  the  perfect  gas  formulation  was  replaced  by 
a  virial  equation  of  state,  generalized  for  a  multi-conqioneat 
mixture  as  summarized  in  Table  V. 

Shock  tube  studies  indicate  the  importance  of  compressibility 
effects  at  higb-pressures.  Figure  5  exhibits  a  10/1  shock  tube 
solution  for  three  different  pressure  ratios  (10/1  atm,  100/10  atm 
and  1(X)0/100  atm).  Note  die  differences  in  the  diock  propaga¬ 
tion  pattern  at  the  higher  pressures.  Shock  speeds  are  fii^  at 
higher  pressures.  The  st^y-state  case  of  the  previous  subsec¬ 
tion  was  studied  again  with  the  new  RAM  version  of  CRAFT. 
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Table  IV.  Dynamic  Grid  Computations 


*  / 


EQUATIONS 

/^Qdv  -  /^Qdv  ♦  0^11  •  ?  dsdt  -  /^‘‘/^^♦dvdt 

Term  I  Tenn  H  Tenn  in 

Tenn  I:  PRIMARY  CONSERVED  VARIABLES 

Term  H:  FLUX  TERMS 

Term  UI:  SOURCE  TERMS 

V(l):  TIME  VARYING  VOLUME 

METHODOLOGY 

•  TREAT  GRID  AS  PURELY  GEOMETRIC  QUANTITY,  i.e..  DEFINITION  OF  CONTRA- 
VARIANT  VELOOTY  UNCHANGED 

•  EVALUATE  FLUX  TERM  BY  ASSUMING  GRID  IS  HELD  CONSTANT  AT  TIME-AVER- 
AGED  VALUE 

•  INCLUDE  ADDmONAL  FLUX  ARISING  FROM  CELL  FACE  MOTION 

•  EVALUATE  SOURCE  TERM  FOR  THE  GRID  AT  THE  NEW  TIME  LEVEL,  THUS 

GIVING: 

-  Q*|y..|  ^  -  V*J  ^ 

F,  ■  fl#,  end  S  is  the  time  avenged  metric 

•  SPECIFY  GRID  MOVEMENT  AS  A  FUNCTION  OF  TIME 

«  EVALUATE  F.  AS  A  THIRD-ORDER  UPWIND  BIASED  FLUX 

•  SOLVE  UNEARIZED IMPUCIT  OPERATOR  USING  ADI  PROCEDURE 

•  USE  NEWTON  ITERATION  FOR  HIGHER  ORDER  ACCURACY  TIME 

Table  V.  Description  of  Gas  Upgrades 


THE  VDUAL  EQUATION  OF  STATE  WAS  INCORPORATED  FOR  A  MULTI-COMPONENT  GAS  MIXTURE 

p-Rr£.2ifi*B.(D£  p,*c.(t)E  (p«n 
•  ’•Vl  ■  •  J 

THE  SECOND  AND  THIRD  VIRIAL  COEFnCIENTS  FOR  THE  MIXTURE  ARE  RELATED  TO  THE  INDIVIDUAL 
SPECIES  VIRIAL  COEFFiaENTS. 

IN  GENERAL 


HOWEVER.  IF  DATA  IS  UNAVAILABLE  FOR  CROSS  COEFFICIENTS  SUCH  AS  B,j,  AN  APPROXIMATE 
RELATIONSHIP  IS  ASSUMED. 


WHERE  Y,  IS  THE  MASS  FRACTION 


Figures  6  and  7  show  the  computed  contours  of  pressure  and 
Mach  number  —  the  shocks  ate  stronger  and  more  closely 
q>aced.  Figures  8  and  9  compare  the  predictions  of  pressure 
along  the  tube  wall  and  projectile  afterbody  from  the  two 
computations.  The  effect  of  compressibility  is  evident  and 
.would  be  more  pronounced  if  combustion  were  included. 

Adaptive  Grid  Considerations 

The  accurate  prediction  of  transient  ram  accelerator  flow- 
fields  is  strongly  keyed  to  the  ability  to  dyi»mically  adapt  the 
grids  to  resolve  the  embedded  shock  discontinuities,  wall 
boundary  layers  and  separated  shear  layer  regions.  Smart  logic 
is  required  to  have  grids  evolve  with  the  transient  structure. 
Research  has  been  initiated  in  this  direction  using  solution 
adaptive  gridding  based  on  the  NASA  Anses  SAGE  adaptive 
solver.^  Applications  of  adaptive  gridding  with  tran- 
sient/unsteady  NS  computations  require  significant  experience 
with  regard  to: 

a)  when  to  implement  the  adaptive  solver,  e.g., 
after  bow  many  iteration  steps; 

b)  the  associated  adaption  parameter,  e.g.,  a 
certain  weighted  combination  of  pressure, 

Mach  number,  temperature,  etc.;  and, 

c)  the  control  parameters,  e.g.,  the  degree  of  grid 
distortion,  orthogonality,  etc.,  per  adaptive 
pass,  etc. 

The  development  of  a  suitable  grid  adaption  strategy  for  ram 
accelerators  bu  been  initiated  with  focus  directed  towards 
developing  rules  for  some  of  the  considerations  outlined  above. 
Pteliminary  adaption  strategies  are  being  developed  for  a  steady- 
state  trajectory  point  analysis  and  will  be  sub^uently  imple¬ 
mented  for  transient  calculations  in  a  *bands-off*  fashion. 
Figures  10  and  11  display  examples  of  adapted  grids  generated 
for  the  steady-state  case  analyzed  earlier.  The  grids  were 
generated  by  selecting  varied  adaption  parameters  and  controlling 
the  grid  distortion.  Figure  12  shows  the  grid  and  resultant 
solution  after  several  passes  of  the  solution  adaptive  grid 
generator  and  further  solution  convergence  with  the  updated 
grids.  Compare  the  crispness  of  the  shock  waves  with  the 
baseline  grid  predictions  of  Figs.  3  and  4. 

Finite-Rate  Chemistry  Considerations 

Finite-rate  chemistry  mechanisms  and  rates  are  required  for 
CFD  simulations.  Unfortunately,  the  vast  majority  of  hydrocar¬ 
bon  kinetics  data^^*^^  were  acquired  (and  mechanisms  formulat¬ 
ed)  under  laboratory  conditions  with  limited  (if  any)  applicability 
,  to  the  high  pressure,  semi-detonative  ram  accelerator  environ¬ 
ment.  Extrapolation  of  kinetics  data  to  higher  pressures  may 
lead  to  numerical  predictions  of  dubious  quality.  Under  such 
constraints,  global  multi-step  mechanisms  provide  a  useful 
compromise  until  more  elaborate  mechanisms  can  be  formulated 
via  the  use  of  computational  chemistry,  as  was  demonstrated 
recently  for  NASP  m  the  hypersonic  scramjet  environment  by 
Cooper  and  coworfcers  at  NASA  Ames. 


One-dimensional  shock  tube  studies  at  laboratory  conditions 
have  been  initiated  to  gain  a  preliminary  understanding  of  the 
kinetic  processes  relevant  to  the  ram  accelerator.  The  objective 
of  these  simulations  is  to  identify  the  interaction  between  shock 
formation  and  ignition  delay  in  a  transient  environment.  A 
related  objective  is  to  also  establish  an  upper  pressure  bound  for 
the  applicability  of  conventional  hydrocarbon  kinetics  rates  and 
mechanisms.  Figure  13  shows  the  schematic  of  a  problem 
studied  and  the  kinetic  mechanism  utilized  for  the  simulations. 
Figures  14  and  IS  show  the  temperature  and  pressure  distribu¬ 
tion  in  the  shock  tube  at  two  instances,  sqparated  by  S/ts.  At  the 
earlier  time,  the  shock  front  is  clearly  seen  (characteriTed  by  the 
von  Neumann  spike),  followed  by  a  constant  pressure  induction 
zone  and  subsequently  die  flame  with  rapid  tenqierature  rise. 
The  flame  front  travels  faster  than  the  shock  and  merges  with  it 
a  little  later,  as  shown  by  the  relatively  uniform  post-shock 
temperature.  Figures  16  and  17  show  the  distribution  of  major 
and  minor  species  at  the  two  time  instances  considered  and  a 
consistent  behavior  is  observed.  Further  one-dimensional  studies 
with  upgrades  rates/mechanisms  are  felt  to  be  essential  for 
gaining  an  understanding  of  what  might  occur  under  actual  ram 
conditions. 

CONCLUDING  REMARKS 

This  paper  has  described  the  CRAFT/RAM  conqwiter  code 
developed  for  application  to  ram  accelerator  flowfidds.  The 
code  developed  is  a  specialized  verrion  of  the  CRAFT  research 
code  with  generalized  thermochemistry,  advanced  turbulence 
modeling  and  specialized  dynrunic/solution  adaptive  gridding 
capabilities.  Issues  pertinent  to  developing  a  fundamental 
understanding  of  the  detailed  physical  and  thmnocbemical 
processes  governing  ram  accelerator  concqits  were  exhibited 
with  numerical  simulatious  illustrating  the  code's  crqrabilities  to 
simulate  these  processes.  Future  work  will  include  boundary 
layer  effects  (projectile  and  launcher  tube),  foil  transient 
simulations,  and  exploration  of  thennochemical  issues. 
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Figure  12.  a)  Grid  obtained  after  multiple  adaptioiif;  b)  Preeaun  contoon  on  adapted  grid 
and  grid;  and  c)  Mach  number  coetoun  oo  adq>ted  grid  and  baaeline  grid. 
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Figure  13.  Reactive  diock  tube  acbemadc  and  methane  Irinrlic  mechaniam. 
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